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The phase diagram of the fully frustrated XY model on a honeycomb lattice is shown to incor- 
porate three different ordered phases. In the most unusual of them, a long-range order is related 
not to the dominance of a particular periodic vortex pattern but to the orientation of zero-energy 
domain walls separating domains with different orientations of vortex stripes. The phase transition 
leading to the destruction of this phase can be associated with the appearance of free fractional 
vortices and is of the first order. The stabilization of the two other ordered phases (existing at lower 
temperatures) relies a positive contribution to domain-wall free energy induced by the presence of 
spin waves. This effect has a substantial numerical smallness, in accordance with which these two 
phases can be observed only in the systems of really macroscopic sizes. In physical systems (like 
magnetically frustrated Josephson junction arrays and superconducting wire networks), the presence 
of additional interactions must lead to a better stabilization of the phase with the long-range order 
in terms of vortex pattern and improve the possibilities of its observation. 
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I. INTRODUCTION 

A uniformly frustrated XY model can be defined by 
the Hamiltonian 



H 



in'} 



(la) 



where the interaction of the variables ipj defined on the 
sites j of some regular two-dimensional lattice is de- 
scribed by periodic function 



Vi9) 



-J cosd , 



(lb) 



J > is the coupling constant, and the summation 
is performed over all pairs of nearest neighbors {jj') 
on the lattice. The non-fluctuating (quenched) variables 
Ajji = —Ajij defined on lattice bonds have to satisfy the 
constraint 



Ajy = 27r/ (mod 27r) 



(Ic) 
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on all lattice plaquettes. The notation □ a below the 
sign of summation implies the directed sum of variables 
Ajj> = —Ajij over the perimeter of plaquette a in the 
counterclockwise direction. For / = the model defined 
by Eqs. ([T|) is isomorphic to the conventional XY model 
(without frustration), whereas for f ~ \ (the maximal 
irreducible value of /, see Ref.'TJ) this model is called fully 
frustrated. 

Since eighties, the uniformly frustrated and, especially, 
the fully frustrated XY models on various lattices have 
been studied rather intensively (for reviews, see Refs. i 
andQ), mostly in relation with experiments on Joseph- 
son junction arrays.^ In these artificial superconducting 
systems variables (pj can be associated with the phases of 
the superconducting order parameter on different super- 
conducting grains forming an array, and Aij is related to 



the vector potential of a uniform magnetic field, whose 
magnitude corresponds to having / superconducting fiux 
quanta per lattice plaquette. The form of Eq. (jlc|) corre- 
sponds to taking into account only the external magnetic 
field and neglecting the field of weak currents fiowing in 
the junctions. Planar magnets with odd number of an- 
tiferromagnetic bonds per plaquette^ are also described 
by the fully frustrated XY models. The recent renewal 
of the interest to the fully frustrated Josephson junction 
arrays has been related to their possible application for 
the creation of topologically protected quantum bits.^ 

In a typical situation, the energy of each of the Joseph- 
son junctions forming an array as a function of the gauge- 
invariant phase difference 



A. 



(2) 



is indeed rather accurately described by the function 
([Tb)) . On the other hand, a magnetically frustrated net- 
work of thin superconducting wires can be described^ 
in the London limit (when the amplitude of the super- 
conducting order parameter can be assumed to be con- 
stant along the wires) by the same Hamiltonian (jlap 
with V{9) replaced by the so-called Berezinskii- Villain 
interaction^'^ having the same symmetry and periodic- 
ity as V{9), see Appendix [A] The fully frustrated XY 
model with the Berezinskii- Villain interaction is more 
convenient for a theoretical analysis, because its parti- 
tion function can be rigorously transformed into that of 
a Coulomb gas with half-integer charges^ (Appendix |A] 
gives more details). In a number of cases it can be ex- 
pected that the main features of the frustrated XY mod- 
els with the conventional and with the Berezinskii- Villain 
interaction are the same. 

The ground states of the fully frustrated XY models 
are characterized by the coexistence of the continuous 
U{1) degeneracy (related to the possibility of the simul- 
taneous rotation of all phases) with a discrete one, whose 
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nature depends on the structure of the lattice. Accord- 
ingly, the fully frustrated XY models allow in addition to 
the phase transition related to unbinding of vortex pairs 
also for the existence of other phase transitions related 
to the disordering of the discrete degrees of freedom. The 
main question in such a situation is what is the sequence 
and the nature of these phase transitions. 

The most thoroughly studied examples of the fully 
frustrated XY models are the models on square and trian- 
gular lattices. In both of them the ground states are char- 
acterized by a regular alternation of the plaquettes with 
positive and negative vorticities. Accordingly, the dis- 
crete degeneracy of the ground states is two-fold, that is, 
the simplest one which is possible. However, it took about 
two decades before it was firmly understood'iiiS why the 
Berezinskii-Kosterlitz-Thouless transition related to un- 
binding of vortex pairs has to take place in these mod- 
els at lower temperatures than the second phase transi- 
tion related to the proliferation of the Ising-type domain 
walls. Due to a strong mutual influence of the two tran- 
sitions (which are situated rather close to each other), a 
convincing numerical demonstration of the Ising nature 
of one of them has required very substantial effortsi^iii 
Before that quite a number of works appeared claim- 
ing that this transition has a non-Ising critical behav- 
ior, which later has been demonstrated to be a finite-size 
effect. 

The properties of the fully frustrated XY model on 
a a honeycomb lattice are not so well understood. The 
family of its ground state is known to have an infinite dis- 
crete degeneracyr^ which can be conveniently described 
in terms of the formation of zero-energy domain walls 
parallel to each other The analysis of the half-integer 
Coulomb gas on a triangular lattice (the dual represen- 
tation of the considered XY model) has revealed^^ that 
in this system the formation of a domain with a different 
charge pattern can cost only a finite energy even when 
the size of this domain is arbitrary large. This has led the 
authors of Ref. [l5| to the conclusion on the absence of a 
long-range order at any nonzero temperature T. 

This conclusion is not directly applicable to the fully 
frustrated XY model with interaction (|lb|) because at 
T > the zero-energy domain walls acquire a positive 
free energy /dw originating from the difference in the 
free energy of the continuous fluctuations (spin waves) 
However, the effect is very weak and according to the 
estimates in Ref. [13 can manifest itself only if the size 
of the system L substantially exceeds Lc ^ 10^, which 
makes it unobservable in real and numerical experiments 
on systems with L Lc. 

Nonetheless, up to now it has remained unclarifled if 
in situations when /dw is absent (as it happens in the 
case of the Berezinskii- Villain interaction) or can be ne- 
glected, the system can still demonstrate a phase transi- 
tion related to the continuous degrees of freedom or the 
destruction of the long-range order of in terms of vor- 
tex pattern prevents the possibility of such a transition. 
The available numerical dat a^^'^^~ — does not lead to a 



unique conclusion and allows for different interpretations, 
including the existence of a spin-glass transitioufi? 

In this work we reexamine the fully frustrated XY 
model on a honeycomb lattice with the aim of finding 
an answer to this question, as well as establishing what 
is the nature of the phase transition (transitions) induced 
by the positiveness of /dw- Our main conclusions can be 
formulated as follows. 

Even when /dw is equal to zero (or can be neglected) , 
this does not mean the complete absence of the long- 
range order at any nonzero temperature. In a such a sit- 
uation, the long-range order is related not to the domi- 
nance of a particular vortex pattern, but to a preferable 
orientation of the zero-energy domain walls. 

The phase transition leading to the destruction of this 
long-range order is related to the appearance of a finite 
concentration of unpaired fractional vortices with topo- 
logical charges q — zL^. The fractional vortices are lo- 
cated at the nodes of the domain-wall network existing 
in the system with /dw = at any temperature, how- 
ever at low temperatures they have to be bound in pairs, 
which imposes the existence of a preferable orientations 
of domain walls. 

In the model with the conventional interaction 
V{9) = — J cos 9 which at T > has small but positive 
/dw (induced by spin- waves), the increase of temper- 
ature leads to the sequence of three phase transitions. 
However, the two of them require for their observation 
very large systems with L ^ Lc ^ 10^. In physical sys- 
tems like magnetically frustrated Josephson junction ar- 
rays and superconducting wire networks, the value of Lc 
will be lower due to the presence of other mechanisms for 
the removal of the accidental degeneracy. 

Although the lowest-temperature phase transition 
(whose critical temperature tends to zero when /dw ~^ 0) 
consists in the appearance of domain walls, it does not 
lead to the destruction of the long-range order in terms 
of vortex pattern. On the contrary, it is related to a par- 
tial restoration of the continuous U{1) symmetry. This 
phase transition belongs to the Ising universality class, 
whereas the two other transitions mentioned above must 
be of the first order. 

A more detailed summary of the phase diagram of the 
fully frustrated XY model on a honeycomb lattice and of 
the properties of different phases and phase transitions is 
presented in the concluding Sec. IVIIII It also compares 
our results with that of other works and discusses their 
relevance for some experimental situations and other uni- 
formly frustrated XY models that allow for the formation 
of parallel zero-energy domain walls. 



II. ZERO TEMPERATURE ANALYSIS 

A. The ground states 

When speaking about both global and local minima 
of the Hamiltonian of an XY model it is convenient to 
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characterize them in terms of vorticities of different pla- 
quettes. Vorticity of plaquette a can be defined as the 
directed sum over the perimeter of this plaquette of the 
gauge-invariant phase differences Ojji = Lpy — ipj — Ajji 
reduced to the interval (— 7r,7r). Here and below we use 
Greek indices for labelling the plaquettes, each of which 
can be associated with a particular site of the triangular 
lattice T dual to the original honeycomb lattice H. 

In the fully frustrated XY model on a honeycomb lat- 
tice variables can acquires values ±7r, ±37r, ±5tt. Since 
the main role in the thermodynamics is played by the 
minima of Hamiltonian ([T]) with Va ~ ±7r (as well as 
by fluctuations in the vicinities of these minima), one 
often speaks not about vorticities but about chiralities 
Ca = Wq/tt = ±1 of different plaquettes. The knowledge 
of the chiralities of all plaquettes is sufficient for restoring 
(up to a global rotation) the values of phase variables ipj 
in the minimum of the Hamiltonian which can be associ- 
ated with this particular configuration of chiralities. 

In all ground states of the fully frustrated XY model 
on a honeycomb lattice the plaquettes with different signs 
of chirality are separated by the bonds with 9jji = ±7r/4, 
whereas the plaquettes with the same sign of chirality are 
separated by the bonds with 9jj' = 0. In Fig. [1] showing 
some of the ground states the plaquettes with positive 
(negative) chirality are shaded (unshaded), whereas the 
bonds with 9jj' = ±7r/4 are shown by arrows. On any 
bond with an arrow, 6jji is equal to +7r/4 if the arrow is 
directed from site j to site j' and to —tt/4 in the opposite 
case. 

Fig. [BJa) shows an example of the ground state with 
the simplest possible structure. In this state the plaque- 
ttes with the same sign of vorticity form parallel straight 
stripes. Accordingly, such states are called below the 
striped states. The set of striped states is characterized 
by a sixfold discrete degeneracy related to the possibili- 



ties of shifting and rotating vortex pattern. Naturally, it 
also has the continuous U{1) degeneracy related to the 
simultaneous rotation of all phases, but in the following 
the term "striped state" will always refer just to the con- 
figuration of chiralities, without specifying the phases. 

Striped states allow for the formation of domain walls 
(separating different realizations of such states) which 
cost no energy. ^'^ An example of such a zero-energy do- 
main wall (ZEDW) is shown in Fig. [Hb). An arbitrary 
number of ZEDWs separated by arbitrary distances can 
be introduced into the system in parallel to each other, 
as is shown in Fig.[ljc). 

The only possible intersection of ZEDWs which does 
not cost extra energy is shown in Fig. [Ud). Since each 
ZEDW can be ascribed a particular direction in accor- 
dance with the orientation of the angle made by stripes, 
as shown in Fig. [1] it is clear that in a ground state the 
system cannot contain more than one intersection of such 
a kind and in the case of the periodic boundary condi- 
tions even one is impossible. Therefore, a typical ground 
state contains an irregular sequence of straight ZEDWs 
which are parallel to each other, as shown in Fig. [IJc). 
The presence of a boundary between two ground states 
with a finite concentration of ZEDWs and different orien- 
tations of these walls would cost an energy proportional 
to the length of this boundary. 

The structure of the ground states is not sensitive to a 
particular form of the interaction in the Hamiltonian of 
the model, as soon as it is even in 9 and behaves more 
or less in the same way as V{9) = —Jcos9 [that is, has 
the single maximum at 6 = tt and the single minimum 
at 6 — whose width is comparable with the period 
of V^(6')]. For example, it can be also the Berezinskii- 
Villain interaction defined by Eq. (IA2p or the piecewise 
parabolic interaction (the zero-tcmperature limit of the 
Berezinskii- Villain interaction) . 




FIG. 1: (Color online) The ground states of the FFXY model on a honeycomb lattice; (a) a striped state; (b) a state with a 
single zero-energy domain wall separating two different striped states; (c) a typical ground state with an irregular sequence of 
parallel domain walls; (d) the only possible intersection of domain walls which does not cost extra energy. The plaquettes with 
positive (negative) chirality are shaded (unshaded). 
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B. Zero-energy domain walls and long-range order 

If domain walls separating different striped states 
would have a positive energy per unit length, at zero 
temperature the system would be frozen in one of the six 
striped states. In terms of chiralities ctq., the order pa- 
rameter corresponding to the corresponding long-range 
order (LRO) can be written as 



a a exp ■ 



(3) 



Here and below denotes the positions of the sites a 
of the dual triangular lattice T which can be associated 
with the plaquettes of the original honeycomb lattice 
whereas Q„ (with 6=1,2,3 and Qi -f Q2 -f Q3 = 0) are 
the three reciprocal vectors of T. In any of the striped 
states one of the three sums Sn = ■5", a grows with 
size of the system as Sn = ±^ {N being the total number 
of the plaquettes in the system), whereas the two other 
are equal to zero. The correlation function 



(4) 



where the average is taken over the six striped vortex 
patterns, has an oscillating behavior. Namely, if triangu- 
lar lattice T is partitioned into four equivalent triangular 
sublattices, then Cchi^a ^ r^) is equal to 1 when a and 
P belong to the same sublattice and to —1/3 when they 
belong to different sublattices. 

When the domain walls shown in Fig. [T] cost no 
energy, the oscillating behavior (without a decay at 
|ra — r^l — > 00) of the correlation function @ (where the 
average now has to be taken over the infinite set of the 
ground states) is retained only when a and j3 belong to 
the same column of sites on 7", that is, when vector Yq^ — vp 
is directed along e„, one of the three lattice vectors of 
T (ei -1-62-1-63 = 0). In such a case, in one third of 
the ground states the line connecting a and /3 cannot be 
crossed by any domain walls and therefore the product 
CTqCt^ depends only on the distance between a and /3. For 
all other directions of — r^, the correlations of CTq, and 
ap are absent, which leads to 



(|5p)cxAr. 



(5) 



This means that the long-range order in terms of Sn,a is 
destroyed. It is worthwhile to emphasize explicitly that 
the dependence ([5]) is not a consequence of algebraic cor- 
relations of (Tq as it could seem from its form, but has a 
different origin explained above. 

However, it is not hard to note that the destruction of 
the long-range order is not complete. In the presence of a 
random sequence of straight parallel ZEDWs, the system 
contains the domains of only four different striped states 
out of six, whereas the two other vortex patterns (with 
stripes parallel to the direction of the domain walls) are 
not represented at all. Therefore, in a typical ground state 
the system still possesses a long-range order. The most 



evident interpretation of this order consists in associating 
it with the direction of ZEDWs. 

Instead of really monitoring the direction of domain 
walls (whose identification requires analyzing chirality 
patterns in rhombic clusters formed by four neighbor- 
ing plaquettes) an order parameter sensitive to the direc- 
tion of domain walls can be introduced by studying the 
distribution of energy between the bonds with different 
orientations. In a typical ground state the average en- 
ergy of the bonds whose direction is perpendicular to the 
orientation of domain walls is equal to F(7r/4), whereas 
for the bonds with two other orientations it is equal to 
i[y(0) + F(7r/4)] < y(7r/4). This allows one to introduce 
an order parameter as 



D 



^exp 



27rn 



jj' ) 



(6) 



where the value of rijji — 1,2,3 depends on the orienta- 
tion of the bond {jj')- 

Thus, the vanishing of the energy of domain walls does 
not lead to the complete destruction of the long-range or- 
der in the system and one can expect that there should 
occur a phase transition related to the destruction of the 
long-range order described by order parameter D. Note 
that D is proportional to the area of the system also when 
one of the six striped vortex pattern is a dominant one. 
However, in contrast to 5, order parameter D does not 
allow to distinguish between two vortex patterns which 
differ from each other by changing the signs of all chiral- 
ities. 

C. Phase correlations 

When speaking about phase correlations, one should 
not forget that a change of gauge leads to the multi- 
plication of exp[i((/3Q, — ipp)] by some phase factor. In 
accordance with that, a gauge-invariant description of 
phase correlations can be achieved by paying attention 
only to the absolute values of the standard correlation 
functions. For taking into account the peculiarities of the 
considered model, it is convenient to introduce the set of 
gauge-invariant phase correlation functions defined as 



Cp{rj - rfc) = \ {exp[ip{ipj - <fik)])\ 



(7) 



and numbered by positive integer p. Another approach to 
introducing gauge-invariant phase correlation functions 
involves considering two identical but completely inde- 
pendent replicas of the system.-- 

If domain walls separating different striped states 
would have a positive energy per unit length, at zero 
temperature there would exist a true long-range order in 
terms of phase variables. In particular, if the honeycomb 
lattice is partitioned into 32 equivalent triangular sub- 
lattices, for any two sites j and k belonging to the same 
sublattice, Cp(rj — r^) would be equal to 1 for any inte- 
ger p. On the other hand, when p is a multiple of eight 
(p = 8p' with integer p'), Cp{rj — r^) would be equal to 
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1 for any pair of sites. This is evident already from the 
fact that for Ojj: — 0,±7r/4 the factors expi(86'jj/) are 
always equal to 1. 

For evident reasons, the presence of a random sequence 
of straight domain walls separating different striped 
states cannot change the behavior of Cp(r) with p — 8p', 
which remain identically equal to 1. However, it leads to 
an exponential decay of Cp(r) for p ^ 8p' for all directions 
except those parallel to e^. For these special directions, 
all phase correlation functions with p < 8 have an oscil- 
lating behavior. Like with chiralities, when one calculates 
the averages of the squares of certain Fourier transforms 
of exp(ip(^j) over all grounds states, they diverge in a 
typical ground state proportionally to the area of the 
system, but this should not be taken for the signature of 
the algebraic correlations of exp(ip93j). 

III. REMOVAL OF THE ACCIDENTAL 
DEGENERACY BY SPIN WAVES 

The infinite degeneracy of the ground states which 
manifests itself through the possibility of the formation 
of ZEDWs is of the accidental nature (i.e., it is not im- 
posed by the symmetries of the system). Accordingly, it 
can removed by the addition of some small interactions 
preserving the symmetries of the Hamiltonian. It is also 
removed at a finite temperature when one takes into ac- 
count the free energy of the continuous fluctuations (spin 
waves). This mechanism for the removal of an acciden- 
tal degeneracyiS^ is often referred to as "order from 
disorder" . In systems with a continuous degeneracy it is 
usually sufficient to compare the contributions from the 
harmonic fluctuations.—"— 

In contrast to that, in the fully frustrated XY model on 
a honeycomb lattice the free energy of the harmonic fluc- 
tuations is exactly the same for all ground states of the 
model. This is a consequence of a hidden symmetry ex- 
isting in the Hamiltonian describing such fluctuations.^- 
The difference in the spin- waves free energy appears only 
when one goes beyond the harmonic approximation. The 
analysis of the third- and the fourth-order anharmonic 
terms reveals that the spin-wave free energy is minimal 
for the striped states whereas the presence of ZEDWs 
adds to the free energy of the system a positive term 
roughly proportional to the total length of the wallsJ^ 
This allows one to ascribe to domain walls the free en- 
ergy per unit length 

/dw = IT'^IJ, 7 « 0.7 X 10-4 , (8) 

where by unit length we mean the lattice constant of the 
dual triangular lattice a a . In terms of a a , the lengths of 
all walls are integer. 

The positiveness of /dw suggests that in the thermo- 
dynamic limit the free energy of a domain wall crossing 
the whole system is infinite and therefore at the lowest 
temperatures (at which the positiveness of /dw is not 
killed by other types of fluctuations) such walls must be 



absent. This ensures the existence of the long-range order 
corresponding to the dominance of one of the six equiv- 
alent striped vortex patterns. On the other hand, at any 
nonzero temperature the presence of spin waves leads to 
an algebraic decay of the correlation functions describing 
phase correlations. 

In terms of chiralities CTq, = ±1 (defined on sites a of 
the dual triangular lattice) the structure of the striped 
state is exactly the same as that of the ground state of 
the antiferromagnetic Ising model on a triangular lattice 
with an additional weak antiferromagnetic interaction of 
second neighbors)^ In this model there exist two differ- 
ent scenarios for the disordering the ordered state^^ the 
choice between which depends on the relations between 
the coupling constants. The basic scenario consist in hav- 
ing a single first-order transition leading to the formation 
of an isotropic domain wall network mixing the domains 
of all six striped states. This leads to the complete loss 
of the long-range order. It is interesting that this phase 
transition takes place when the free energy of a domain 
wall (per unit length) is still positive and not when it 
vanishes. 

However, for some relations between the coupling con- 
stants the first-order transition related to the formation 
of a domain wall network can be preceded by a continuous 
phase transition related to vanishing of the free energy 
of double domain walls. This phase transition leads only 
to a partial loss of the long-range order. Namely, above 
it the threefold degeneracy with respect to the directions 
of stripes still persists, but the system consists of alter- 
nating domains of the two striped states with the same 
orientation of the stripesi^ 

The existence of different options suggests that in the 
fully frustrated XY model on a honeycomb lattice, the 
conclusion on how the disordering of the striped state 
takes place also should be based on the comparison of 
different mechanisms for disordering. In particular, it still 
remains to be established what is the nature of the first 
(with the increase in temperature) phase transition in- 
duced by the extreme smallness of the domain wall free 
energy and whether it leads to the complete loss of the 
long-range order, or only to its partial loss which leaves 
a place for another phase transition (or transitions). 

In order to answer these questions we have to under- 
stand what is the most efficient mechanism for the pro- 
liferation of domain walls. And since the decrease of the 
free energy of domain walls is induced by the presence on 
them of point-like defects we have to find the point-like 
defects which play the dominant role in the development 
of the domain- wall fiuctuations. 



IV. POINT-LIKE DEFECTS 
A. Conventional vortices and vortex pairs 

As any other XY model, the fully frustrated XY model 
on a honeycomb lattice allows for the formation of con- 
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ventional vortices with integer topological charges. In 
particular, a vortex with topological charge g = ±1 is 
created when the sign of chirality of some plaquette is 
flipped in comparison with what it would be in a ground 
state. On any contour surrounding the core of a vor- 
tex with topological charge q the deviations of gauge- 
invariant phase differences Ojji = Lpy — Lpj — Ajj> from 
the values they would have in a ground state (0 or ±7r/4) 
sum up to 27rg. 

At low temperatures vortices with integer topological 
charges can be present in the system only in the form of 
small neutral pairs whose concentration at T J has 
to be exponentially low. Below we demonstrate that all 
phase transitions in the considered model take place at 
T <^ J, when the bound pairs of integer vortices are 
exponentially rare and, accordingly, have no influence on 
the properties of the system. For this reason the presence 
of such pairs can be neglected. 

In addition to conventional vortices, in frustrated XY 
models there also can exist fractional vortices localized 
on corners and intersections of domain walls?^ In situa- 
tion when the energy of domain walls connecting different 
fractional vortices is exactly equal to zero, the energies 
of neutral pairs of fractional vortices will be much lower 
than that of pairs of conventional vortices. Accordingly, 
one can then expect the neutral pairs of fractional vor- 
tices to be the most important finite-energy excitations 
at the lowest temperatures. 

B. Fractional vortices 

In any ground state of the fully frustrated XY model on 
a honeycomb lattice each plaquette has exactly two near- 
est neighbors with the same sign of chirality and four with 
the opposite sign. In such a situation the simplest idea 
of having a point-like excitation consists in constructing 
a configuration in which one (and only one) of the pla- 
quettes has a wrong number of neighbors with the same 



sign of chirality. Indeed, such configurations can be con- 
structed, see Fig. where as in other figures below the 
plaquettes with positive (negative) chiralities are marked 
by the presence (absence) of a shading. 

In particular, in Fig.[2l^a) the plaquette denoted by the 
sign plus has three neighbors with the same sign of chi- 
rality (instead of two) , whereas the analogous plaquette 
in Fig.[21[b) one neighbor. In both configurations, in order 
to have the correct number of neighbors with the same 
sign of chirality on all other plaquettes, three ZEDWs 
have to merge at the defect. 

An important property of these defects is that they are 
fractional vortices with topological charges q = ±1/8. In 
any ground state, the correct value of vorticity on each 
plaquette is achieved because on all four bonds separating 
it from the plaquettes with the opposite sign of vorticity 
variables djj' are equal either to -|-7r/4 or to — 7r/4 and 
sum up to ±7r. When the number of neighbors with the 
same sign of vorticity is equal to one or three (instead of 
two), there appears a misfit of ±7r/4 which has to be dis- 
tributed over all six bonds surrounding the corresponding 
plaquette. The same misfit is also present on any closed 
contour surrounding the defect independently of the size 
of this contour, from where it is clear that the deviations 
of variables 9jjr from the values they would have in a 
ground state (0 or ±7r/4) decay with the distance from 
the defect R as 1/R. This leads to a logarithmic diver- 
gence of the energy of a single defect and a logarithmic 
(at large distances) interaction between such defects. In 
accordance with their topological charges q = ±1/8, the 
logarithmical interaction of the fractional vortices in the 
considered problem is 64 times smaller than that of con- 
ventional vortices with topological charges q — ±1. 

The defect shown in Fig.[2][c) (in which the central pla- 
quette has no neighbors with the same sign of vorticity) 
can be considered as an overlap of two defects of the type 
shown in Fig. HKb). This statement applies both to the 
number of ZEDWs which meet at the defect and to its 
vorticity. 




FIG. 2: (Color online) The simplest point-like defects in which one of the plaquettes has a wrong number of neighbors with the 
same sign of chirality. 
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FIG. 3: (Color online) Neutral pairs of fractional vortices 
which cannot be present in the system in the absence of free 
fractional vortices. 

In Fig.[2]the positions of fractional vortices are denoted 
by the sign plus, because in all configurations shown in 
this figure the topological charge of the fractional vor- 
tex is positive. Below the positions of fractional vortices 
are denoted by pluses and minuses in accordance with 
the signs of their topological charges. In the case of the 
Berezinskii- Villain interaction, it is possible to exactly 
express the energy of any configuration of chiralities in 
terms of the pairwise interaction between fractional vor- 
tices (see Appendix [B|) . 

C. Neutral pairs of fractional vortices 

Now we know how the fractional vortices look like and 
the next step is to find how two fractional vortices can 
be combined to form a neutral pair with a finite energy. 
Since each ZEDW can be ascribed a particular direction 
as shown in Figs. [T]and[5J each fractional vortex can be 
considered either as a "sink" [the configuration in Fig. 
[5]Ja)] or as a "source" [the configurations in Figs, ^h) 
and[5Jc)] of domain walls. 

Apparently in the case of periodic boundary conditions 



the number of sources must be equal to the number of 
sinks [naturally, for the correct balance the configurations 
of the type shown in[21^c) have to be counted as double 
sources]. It is impossible to construct an isolated pair of 
fractional vortices from two sink configurations because 
the domain walls meeting in them would have to inter- 
sect each other, which would imply the existence of other 
defects close by. For the same reason it is impossible to 
construct an isolated pair from two double sources. On 
the other hand, although locally a neutral pair consist- 
ing of two sources looks like a legitimate object (see Fig. 
[3]), such pairs cannot appear because in the absence of 
unpaired fractional vortices there will be no sinks to com- 
pensate for these sources [in that respect they resemble 
the configuration shown in Fig. [IJd)]. 

Therefore, in a domain wall network formed by small 
pairs of fractional vortices, each pair should consist of a 
source and a sink. Such pairs can be divided into three 
classes, exemplified in Fig. The first of them is an in- 
tersection of a single ZEDW and a double domain wall, 
that is, a pair of parallel ZEDWs (which change their 
orientation after crossing a single wall). An example of 
such a pair of fractional vortices is shown in Fig. |31[a). 

The second type of the fractional- vortex pairs are the 
bends on a double domain wall, see Fig. |H[b) . For such 
a pair to be neutral the distance between the two sin- 
gle walls forming the double wall has to be even (in the 
considered model it is natural to measure the distance 
between domain walls in units of /ia, the height of the 
triangular cell of the dual lattice). If the distance between 
two parallel ZEDWs would be odd (for example, the min- 
imal distance), in the configuration analogous to the one 
shown in Fig. Hl^b) both fractional vortices would be of 
the same sign. 

The third class of the fractional- vortex pairs corre- 
sponds to a double domain wall ending on a single domain 
wall, see Fig. Ufc). In this case the distance between the 
two parallel ZEDWs also has to be even in order to avoid 
a situation when both fractional vortices in a pair are of 
the same sign. 
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V. THE LOWEST-TEMPERATURE PHASE 
TRANSITION 

When a linear defect (for example, a domain wall) has 
proper energy or free energy per unit length e and addi- 
tionally allows for the creation of point-like defects with 
energy Ed, its total free energy per unit length is given 
by 

F ^ e - cT\n (^1 + e-^""/^^ K e ~ cTe''^^''^ , (9) 

where c is the number of places per unit length available 
for the formation of a point-like defect and we have as- 
sumed that T <C -Ed- A natural choice for the unit length 
when speaking about domain walls in a lattice model is 
the lattice constant either of the direct or of the dual 
lattice, depending on whether the domain walls can be 
associated with the bonds of the former or of the latter. 
Then c has to be equal either to one or to some simple 
fraction (or integer) of the order of one. In our system, 
it is convenient to measure the length of domain walls in 
lattice constants of the dual (triangular) lattice, ca- 

It follows from Eq. (jQ]) that for e <C i?D the temper- 
ature at which F changes sign lies between e/c and 
(e/c ^ T:,, <C E-o) but depends on E-q much stronger 
than on e: 



Et 



ln(cr*/e) 



(10) 



It is clear from Eq. ([T0| that the first linear defects whose 
free energy vanishes with the increase in temperature are 
those that allow for the formation of point-like defects 
with the lowest energy. 

It follows from the analysis of Sec. lIV Cl that the lowest- 
energy point-like excitation on a linear defect in the con- 
sidered problem is a bend on a doubly-spaced double do- 
main wall (DDDW), consisting of two single domain walls 
separated by distance 2/ia, see Fig. HJ^b). A numerical 
minimization of energy in a finite system (complemented 
with the extrapolation to the infinite size) shows that in 
the model with V{9) = — J cos 6 the energy of this defect 
Eb is one order of magnitude smaller than the coupling 
constant, Eb/J « 0.111. Taking for DDDW e ^ 2/dw 



with /dw given by Eq. ([8]), c = 1 and i?D = Eb, one ob- 
tains from Eq. (fTO|) that TdddWj the temperature where 
the free energy of a DDDW vanishes, is about 14 times 
lower than Eb and two orders of magnitude lower than 
the coupling constant, 



Tdddw ~ 0.81 X 10"^ J. 



(11) 



In the vicinity of TdddWi the influence of the point- 
like defects on the free energy of other types of domain 
walls can be neglected. This follows from the analysis of 
their energies. In particular, a bend on a double domain 
wall with the minimal separation between the walls, /ia 
(a singly-spaced double domain wall), contains two frac- 
tional vortices with the same sign of topological charge 
q = ±i. Therefore, its energy is logarithmically diver- 
gent and the lowest-energy excitation on such a wall is 
not a bend but a kink formed by two bends, see Fig.lSJa). 
Since this object can be considered as a neutral pair of 
fractional vortices with topological charges q = ±|, its 
energy E^ has to be larger than that of a bend on a 
DDDW (a neutral pair of fractional vortices with topo- 
logical charges q = ±|)- This is confirmed by a numerical 
calculation revealing that Ek ~ 0.287 J, that is, the en- 
ergy of such a kink is 2.6 times larger than that of a 
bend on a DDDW. This ensures that at T w Tdddw the 
correction to the free energy of a singly-spaced double 
domain wall is nine orders of magnitude smaller than its 
bare value and therefore can be neglected. 

A single domain wall even cannot make a bend, 
because the direction of such a wall is uniquely de- 
termined by the directions of vortex stripes in the 
states which it separates. The simplest point-like 
excitation on a single domain wall is a kink with 
height 2/iA- Such a kink is formed by four fractional 
vortices with q — ±1/8 having the same sign [see Fig. 
[5{b)] and therefore has topological charge q = ±i. 
Accordingly, the minimal- energy excitation on a single 
domain wall is a neutral pair of such kinks forming 
a kink of height 4ft,A, see Fig. [DJc). It is clear that 
the energy of this complex object has to be about 
one order of magnitude larger than that of a bend on 
a DDDW (for the case of the Berezinskii- Villain interac- 




FIG. 5: (Color online) Kinks on a singly-spaced double domain wall (a) and a single domain wall [(b) and (c)]. 
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tion, the values of Eb, Ek, E^^ and other finite-energy 
defects mentioned in the text can be found in Appendix 

El. 

From this one can conclude that in the vicinity of 
Tdddw all other topological excitations save DDDWs can 
be neglected. Note that since the vortex patterns on the 
two sides of a DDDW are the same, a DDDW is not a 
topologically stable object and in principle can end some- 
where, see Fig. inija). This is especially evident when one 
notices that any DDDW is nothing else but a line of pla- 
quettes with alternating chiralities on which the signs of 
all chiralities are reversed. However, the states on the two 
sides of a DDDW are not exactly the same, namely, they 
differ by a phase rotation by tt. In order to compensate 
for this misfit an end-point of a DDDW has to be a frac- 
tional vortex with topological charge q = ±i. The same 
conclusion also follows from observing that an end-point 
is formed by four fractional vortices with g = ±i having 
the same sign, see Fig. (Blja). 

At temperatures we are discussing, vortices with topo- 
logical charges q — ±^ have to be bound in small neutral 
pairs. Therefore, in the nearest vicinity of any end-point 
there always has to be present another end-point with the 
opposite topological charge, so that the two DDDWs can 
be considered as a continuation of one another [analogous 
situation is known to exist in the unfrustrated XY model 
in which the interaction function V{9) in addition to the 
main minimum at ^? = has an additional minimum at 
6 ~ TT with almost the same dcpth^^]. The orientation 
of a fluctuating DDDW is determined by the orientation 
of stripes in the vortex pattern on its sides, on the aver- 
age such a wall will be perpendicular to the direction of 
stripes. 

At T > Tdddw the free energy of a single DDDW be- 
comes negative which suggests that there should appear 
a finite concentration of such walls p restricted by their 
repulsion. Since this repulsion is of contact nature, the 
dependence of p on the distance from the critical temper- 
ature can be expected to be of the Pokrovsky-Talapo\s2^ 
type, p(x {T - Tdddw)^/^- 

However, in addition to the formation of infinite 
DDDWs the system allows also for the creation and an- 
nihilation of pairs of fluctuating DDDWs, see Fig. (Bl^b). 
This means that in addition to infinite DDDWs crossing 
the whole system, there should be present finite-size de- 
fects formed by two fluctuating DDDWs, which on both 
sides of the defect merge together as shown in Fig. (Hb). 
In the theory of the commensurate- incommensurate tran- 
sitions the objects where n solitons merge together are 
known as dislocations and the case of n = 2 corresponds 
to the double degeneracy of the ground state.-^ In our 
system, the presence of a two-fold degeneracy manifests 
itself through the fact that the states on the two sides 
of a DDDW differ from each other by phase rotation by 
TT. As a consequence of this, after crossing two DDDWs 
one returns to the same state as before. Accordingly, the 
phase transition related to the proliferation of DDDWs 
has to be of the Ising type. 



Since the states on the two sides of a DDDW corre- 
spond to the same vortex pattern but differ by a phase 
rotation by tt, the phase transition related to their prolif- 
eration does not destroy the long-range order in terms of 
chirality but leads to a partial suppression of phase cor- 
relations. Namely, an algebraic decay of correlation func- 
tion Ci(r) is replaced by by an exponential one, whereas 
the correlations of the double phase [described by C2(r)] 
remain algebraic. In the Pokrovsky-Talapov regime, the 
temperature dependence of the correlation radii describ- 
ing the decay of Ci(r) is characterized by two different 
values of the exponent i/ {ly — 1/2 and ly — 1) ioi the two 
directions, whereas in the Ising regime the value of = I 
is the same for both directions. 

The crossover from the Pokrovsky-Talapov behavior 
to the Ising one must take place very close to the phase 
transition temperature, because each point where two 
DDDWs are created or annihilate is a pair of fractional 
vortices with topological charges q = zt^ and accord- 
ingly has larger energy than a bend on DDDW. It should 
be emphasized that when the temperature is about 14 
times smaller than Eb even a relatively small increase 
in energy in comparison with Eb leads to the suppres- 
sion of the Boltzmann factor exp(— i?D/T) by orders of 
magnitude. 

Just above Tdddw the distance between the DDDWs 
L is much much larger then their "effective width" ^ given 
by the average distance between neighboring bends on a 
wall, ^ = £,{T) — exp{EB/T). Since it is evidently disad- 
vantageous to have L ^ ^ (this would force the concen- 
tration of bends with positive energy to be much larger 
than the optimal one), the ratio L/S, has to saturate with 
the increase in temperature at L/^ ~ 1. For this to take 
place, the Boltzmann factor 

wb{T) = exp{-EB/T) (12) 

has to become much larger than fuw/T, which allows one 
to neglect the first term in Eq. ([9]). Due to the exponential 
dependence of wb(T) on T, this is achieved already in a 
close vicinity of Tdddw- 

In this regime the free energy of the fiuctuating 




(a) (b) 

FIG. 6: (Color online) (a) an end-point of a DDDW is a g = 
±i vortex; (b) an end-point of a finite size defect formed by 
two fluctuating DDDW. 
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DDDWs (per lattice plaquette) is of the order of 
— rexp(— 2i?B/T), where one factor exp(— Eb/T") is di- 
rectly related to the energy of the elementary point-like 
defect and the other one appears because the distance be- 
tween neighboring DDDWs has to be exponentially large 
and, accordingly, the number of places where these point- 
like defects can be created is exponentially suppressed. 
However, in the situation when the fluctuation-induced 
free energy of domain walls can be neglected (or is ab- 
sent from the beginning, as in the case of the Berezinskii- 
Villain interaction), one has to consider also other pos- 
sibilities for the appearance of point-like defects in the 
system. This task is carried out in the next section. 

VI. UNIAXIAL NETWORK STATE 

In this section we first discuss what happens at low 
temperatures if /dWi the free energy of straight domain 
walls per unit length, is exactly equal to zero and there- 
fore does not induce a removal of the accidental degener- 
acy of the ground states. And after finding an answer to 
this question, we return to the problem with small but 



finite /dw > 0. 

For /dw — 0; 0116 can expect that at low tempera- 
tures a typical configuration in addition to having paral- 
lel straight domain walls may include a small concentra- 
tion of local defects (with finite energies) at which these 
domain walls change their orientations. However, a finite 
concentration of such objects can appear at arbitrary low 
temperature only if the entropy (per defect) of the net- 
work of domain walls connecting them can be made ar- 
bitrary large by the decrease of their concentration, as in 
the gas of non-interacting particles. Without a detailed 
analysis it is not evident what structure such a network 
can have and whether its appearance induces a complete 
disordering of the vortex pattern. 

In Sec. II VI we have shown that the lowest-energy local 
defects in our system are pairs of fractional vortices with 
topological charges q = ± i . Four different classes of such 
defects are illustrated in Figs. |3] and ID However, a net- 
work of domain walls cannot contain the pairs of the type 
shown in Figs. [3] because the domain walls "emitted" by 
them can end only on free (that is, unpaired) fractional 
vortices, which at low temperatures cannot be present in 
the system. 




FIG. 7: (Color online) A network formed by parallel single domain walls connected by double-spaced double domain walls: 
(a) with relatively small density riDw of single walls ; (b) with now ~ 1/2. 
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The neutral pair of fractional vortices shown in Fig. 
mja) is an intersection of single and double domain walls. 
In order to have a finite concentration of the defects of 
such a kind, it is necessary to have a sequence of paral- 
lel single domain walls and a sequence of double domain 
walls crossing them. However, the entropy of such a net- 
work will be proportional to the number of the walls in 
these sequences, which does not give the entropic contri- 
bution to free energy a chance to overcome the positive 
term related to the proper energy of the defects and pro- 
portional to the number of the intersections. This sug- 
gests that such defects cannot play a substantial role in 
the spontaneous formation of a domain-wall network at 
the lowest temperatures. 

In contrast to that, the neutral fractional- vortex pairs 
of the two other types shown in Fig. Sfb) and Fig. UJ^c) 
allow for having an arbitrary large entropy per defect. In 
particular, the configuration shown in Fig. Sfb) is noth- 
ing else but a bend on a double-spaced double domain 
wall. In the last paragraph of Sec. |V] we have argued 
that for /dw the free energy (per site) of the se- 
quence of such double walls has to be of the order of 
— Texp(— 2_EB/r), where Eb is the energy of a single 
bend. In this estimate, the exponent contains not just 
the ratio E-q/T like it would be if different point-like de- 
fects could be created independently of each other, but 
2E'q/T. The second factor exp(— iJe/T") appears because 
double domain walls have to be separated by distances of 
the order of ^ = exp(i?B/T'), and therefore the number 
of places where the point-like defects with energy E^, can 
appear is exponentially suppressed. 

Another possibility to get a free energy of comparable 
order consists in considering configurations in which the 
energy of point-like defects is larger than i?B, but the 
number of places available for their creation is just pro- 
portional to the area of the system without an additional 
exponential suppression. This happens when one consid- 
ers a configuration with a finite density p\\ of parallel 
single domain walls and inserts into it a finite concentra- 
tion of point-like defects each of which is a double-spaced 
double domain wall starting on one single wall and end- 
ing on a neighboring one, as shown in Fig. [7l[a). The 
energy of these defects is close to 2£'b and, accordingly, 
for p|| <C 1 the free energy of such a "uniaxial" network 
can be estimated as 

i^UNW - r exp(-2SB/T) . (13) 

Note that Fig. [7] gives just a schematic representation 
of the structure of corresponding states. In reality, at 
T <C E-Q the typical distances between neighboring local 
defects have to be exponentially large. Accordingly, the 
mutual influence of the local defects can be neglected. 

It is clear from Eq. (fT3|) that it is more profitable to 
have p\\ of the order of one rather than p\\ <C 1. The 
concentrations close to 1/2 (the average concentration of 
parallel domain walls in a typical ground state) seem to 
be the most optimal because they optimize both the num- 
ber of the available configurations as well as the energy 



of typical local defects. 

The minimum of energy of two fractional-vortex 
pairs connected by a double-spaced double domain wall 
is achieved when the four fractional vortices form a 
symmetric rombus, as shown in Fig. [5]Ja). Both for 
V{6) = — JcosO and for the Berezinskii- Villain interac- 
tion, the energy of such a defect, E^ is smaller than 2Eb 
by approximately 11%. The same is true for another de- 
fect with the same rhombic arrangement of four fractional 
vortices shown in Fig. [8l^b). Accordingly the free energy 
of the domain wall network containing such defects (and 
maybe some other defects as well) at the lowest temper- 
atures will be of the order of —Texp{—Eii/T), that is, 
lower than the free energy of the sequence of fluctuating 
double walls. 

Note that the list of local defects whose energy is lower 
than Eji in addition to configurations which cannot par- 
ticipate in the formation of a domain-wall network [Fig. 
[3] and bends on double-spaced double domain walls [Fig. 
Ufb)] includes only the intersections of a double domain 
wall either with a single wall [Fig. ^a) ] or another dou- 
ble wall [Fig. [9] . The presence of such local defects is not 
sufficient for the construction of a domain wall network 
whose entropy dominates over the energy of the defects. 
This ensures that for /dw = Oi the main role in the 
formation of the domain-wall network is played by the 
rhombic defects [shown in Fig. [8] whose presence leads 
to the formation of the uniaxial network state schemati- 
cally depicted in Fig.[7l^b). 

Reminder that in a typical ground state the concen- 
tration of domain walls parallel to each other is equal to 
1/2. For /dw = at the lowest temperatures the system 
will be in the same phase, the main difference consisting 
in the presence of an exponentially small concentration 
of local defects, which shift the positions of domain walls 
but do not lead to the change of their orientation, see 
Fig. [7{b). Like at T = 0, the presence of a finite concen- 
tration of parallel domain walls leads to the intermixing 
of four out of six striped vortex patterns. However, there 
remains a triple degeneracy related to the orientation of 
domain walls. The order parameter suitable for describ- 
ing this kind of ordering has been introduced in Sec. Ill Bl 




(a) (b) 

FIG. 8: (Color online) The structure of the typical local de- 
fects participating in the formation of the uniaxial network 
state [shown in Fig. [TJb)]. 
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In terms of phase correlations, the main qualitative 
difference with the case of T = is that at T > cor- 
relation function Cp(r) with p — 8 acquires an algebraic 
decay (analogous correlation functions with p < 8 decay 
exponentially even at T = 0). However, this property 
is related not to the appearance of local defects in the 
domain-wall network but to the presence of spin waves. 

Thus, we have found that the formation of local defects 
works as a mechanism for the removal of the accidental 
degeneracy of the ground states which makes the free 
energy of the system dependent on the concentration of 
parallel single domain walls with F having the 
minimum in the vicinity of p|| = 1/2 (in the limit of 
T — > exactly at pn = 1/2). Now we have to remem- 
ber that in the absence of this mechanism i^(p||) is equal 
at small to /dw(T')P|| and, accordingly, is minimal at 
P|| = 0. When both mechanisms are taken into account, 
the function F at least in some interval of temper- 
atures will have two minima whose depths change with 
temperature. At the temperature where these two depths 
become equal, a first-order phase transition has to take 
place. 

For /dw > 0, a possibility to describe the system as- 
suming /dw — exists if the ratio 



/dw(T) 



Texp{~En/T) 



(14) 



is much smaller then one. For the model with 
V{9) = — J cos 9 the substitution of the expression ([8]) for 
hw{T) into Eq. (HH) gives that at T/J = 0.02 the value 
of k{T) is close to 0.03. Accordingly, at T/J « 0.02 the 
corrections to free energy related to the positiveness of 
/dw can be neglected and one can expect that the free 
energy of the system is minimal in the uniaxial network 
state. On the other hand, at T/J ^ 0.015 the negative 
contribution to the free energy related to the appearance 
of the point-like defects with energy E^i cannot overcome 
the positive term related to the proper free energy of do- 
main walls and any reasons for the appearance of a uniax- 
ial network with a large concentration of parallel domain 
walls are absent. 

From this it is clear that the first-order transition men- 
tioned above takes place at some temperature between 





(a) (b) 

FIG. 9: (Color online) Intersections of double domain walls. 



0.015 J and 0.02 J. Below this transition the minimum 
of free energy is achieved in the phase with a sequence 
of fluctuating double walls described in Sec. |V] and above 
it in the uniaxial network with a relatively high density 
of parallel domain walls. The conclusion that this phase 
transition is of the first order is additionally confirmed by 
checking that the presence of fluctuating double-spaced 
double domain walls does not lead to a decrease of free 
energy of a single domain wall crossing them. The ad- 
vantages of the uniaxial domain-wall network manifest 
themselves only when the density of single domain walls 
is of the order of 1/2. They consist in the possibility of 
having local defects of an additional type shown in Fig. 
[51 which increases the number of configurations with the 
given number of the defects, that is with the given energy. 



VII. PHASE TRANSITION TO A DISORDERED 
PHASE 

For the complete disordering of vortex pattern the net- 
work of domain walls has to have a structure leading to 
an unbiased representation of all six striped patterns. If 
all fractional vortices are bound in small neutral pairs 
this is impossible. The existence of a finite concentration 
of such pairs as those shown in Fig. [3] is prohibited by a 
simple conservation law explained in Sec. IIV CI whereas 
the presence of small concentration of pairs of the three 
types represented in Fig. 2] leads to the intermixing of 
only four striped states out of six. Such a uniaxial net- 
work has a preferable orientation of domain walls. 

In order to have an equal representation of all six 
striped vortex patterns there should exist a finite con- 
centration of free fractional vortices which are not bound 
in pairs. If their positions could be arbitrary (that is, not 
restricted by the requirement for them to be connected 
by domain walls) free fractional vortices would appear 
when the two logarithmically divergent contributions to 
the free energy of a single fractional vortex 



FpY — EpY — TSpY 

compensate each other. Here 



Epy Eq In L„ 



Eq 



(15) 



(16) 



is the energy of a fractional vortex and Spy ~ 2 In Lmax 
its entropy, whereas F = F(r) is the helicity modulus^ 
and Lniax ^ oo the size of the system. The prel- 
ogarithmic factor in .Ffv becomes equal to zero at 
T = ^Eq = -^T. On a honeycomb lattice with the con- 
ventional or the Berezinskii- Villain interaction F(r = 0) 
is respectively slightly below or slightly above J/2. Sub- 
stitution of this value suggests that in the absence of any 
restrictions related to domain walls free fractional vor- 
tices would appear at Tpv ~ O.OlJ.-^i 

Naturally, in the fully frustrated XY model on a honey- 
comb lattice the entropy of the system of fractional vor- 
tices is substantially decreased by the requirement that 
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they have to be connected by straight domain walls (with 
maybe some rare kinks on them). Moreover, the direc- 
tion of each domain wall is uniquely determined by the 
directions of the vortex stripes in the two states which 
it separates. Nonetheless, these rather strong restrictions 
still allow for constructing a domain- wall network which 
leads to an unbiased representation of all six striped vor- 
tex patterns. A possible structure of a such network is 
schematically shown in Fig. [TOl Here the letters A, B 
and C are used to denote the domains with three differ- 
ent orientations of stripes. Note that all walls between A 
and B are parallel to each other. The same is true for all 
walls between B and C, as well as for all walls between 
C and A. For the sake of lucidity Fig. [10] does not show 
which of the two versions of A, of B, or of C (related to 
the change of the signs of all chiralities) is realized in each 
particular domain. This depends on the exact positions 
of domain walls. 

Quite remarkably, the entropy (per node) of such a 
network logarithmically depends on typical distance be- 
tween neighboring nodes^ like it would be in a gas of free 
particles. This is so because each domain of a network 
with such a structure can be shifted in parallel to six do- 
main walls which end up in its corners without changing 
the number of domains or the total length of the domain 
walls. Thus, if the typical distance between neighboring 
nodes of the network is of the order of L ^ 1, then each 
domain can occupy a number of positions which is pro- 
portional to L. Since the number of domains is equal to 
one half of the number of nodes, the domain- wall-network 
entropy per fractional vortex is given (for large enough L) 
by 5 In L per fractional vortex. This quantity is only four 
times smaller than in the absence of restrictions induced 
by the presence of domain walls. 

The logarithmic behavior of the domain-wall network 
entropy was first discovered by Villain^ when consider- 
ing a honeycomb network in which each domain has the 
shape of a hexagon with all angles equal to 120 degrees. 
Such a network is formed in a system with a threefold 
degeneracy in which a domain wall of a given type (for 
example, a wall between states A and B) can have only 
three particular orientations out of six that seem to be 
possible when one looks just at the symmetry of the lat- 
tice. A domain- wall network with the structure shown 
in Fig. [10] (in which each domain also has the shape of 
a hexagon but with angles of 60, 120 and 240 degrees), 
has been proposed for the Ising model on a triangular lat- 
tice with the antiferromagnetic interaction of the nearest, 
second and third neighbors. For some relation between 
parameters, the ground states of this Ising model have 
the striped structure^- completely analogous to that of 
the striped ground states of the considered XY model. 
However, the important difference between the two mod- 
els is that in the model we discuss now the nodes of the 
domain wall network have a long-range logarithmical in- 
teraction. 

In accordance with that, in the fully frustrated XY 
model on a honeycomb lattice the free energy (per lattice 



plaquette) of a domain wall network with the structure 
shown in Fig. [10] is given by 

Fnw = ^ (^0 - I) In ^ + F^Ul) , (17) 

where ci '-^ 1 is a constant of purely geometrical origin, 
whereas L ^ 1 is the typical distance between neighbor- 
ing fractional vortices. The last term in Eq. P7|) describes 
the dominant corrections to the logarithmic terms and for 
L ^ 1 can be estimated as F^^{L) ~ Eq/L"^. 

At T < Tfv = 2Eo, when both terms in Eq. ^ are 
positive the minimum of -Fnw is achieved when L oo, 
that is, in the absence of a domain-wall network. On the 
other hand, at T > Tpv the free energy of the system at 
large L decreases with the decrease in L. Therefore, the 
optimal value of L at such temperatures is determined by 
the interplay between the two terms in Eq. (jl7p . It is then 
clear without performing any calculations, that when T 
tends to Tpv from above, the value of L minimizing -Fnw 
goes to infinity. 

However, all the conclusions of the previous paragraph 
would be applicable only if any other fluctuations save the 
formation of a domain wall network with the structure 
shown in Fig. [TU] are prohibited. In reality, one expects 
that at such temperatures the system is in the uniaxial- 
network state with a large concentration of parallel do- 
main walls. In this situation, the structure of the domain- 
wall network related to the appearance of unpaired frac- 
tional vortices also can be illustrated by Fig. [TO] where 
now the domains denoted by the same letter correspond 
to uniaxial-network states with the same dominant ori- 
entation of domain walls. Nonetheless, exactly as before, 
each node of the network corresponds to an unpaired 
fractional vortex with topological charge g = ±i, which 
means that one can still use expression (|17p for describing 
the energy of the interaction of unpaired fractional vor- 
tices plus the entropy of the network related to fluctua- 
tions of positions of different domains. However, one also 
has to add to this expression the energy of the boundaries 
between different uniaxial networks. 

A boundary between two uniaxial-network states 
with different preferable orientations of domain walls is 




FIG. 10: The structure of a domain wall network in which all 
six regular vortex pattern are intermixed. 
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schematically shown in Fig. I VIII where each line corre- 
sponds to a domain wall, whereas letters A, B and C 
again denote domains with three different orientation of 
stripes. Note that one of the uniaxial networks (the one to 
the right) contains mainly domains A and C, whereas the 
other one domains B and C. It is clear that the change 
of the dominant orientation of domain walls which has 
to take place at a boundary between domains requires to 
have on this boundary a large concentration of fractional 
vortices and accordingly the energy of such a boundary 
per unit length has to be of the order of Eq. 

Accordingly, the expression for the free energy of the 
domain wall network formed by unpaired fractional vor- 
tices separated by distances of the order of L, must also 
contain a positive contribution of the order of Eq/L, 
which for large L will always dominate over the first term 
in Eq. (ITTl) . This suggests that the phase transition re- 
lated to the appearance of free fractional vortices must 
take place at T = Tdis > Ipy and has to be of the first 
order, because the two minima of the function F^^{L) 
(at L = oo and at finite L) will always be separated by 
some barrier. Moreover, since the behavior of the function 
Fnw{L) is determined (in addition to T) by the single 
energy scale, Eq, one can expect that the shift of the 
transition temperature, Tdis, with respect to Tpv will be 
of the order of Eq, whereas the typical distance between 
fractional vortices just above the transition can be only 
numerically larger than one, and not parametrically. Ac- 
cordingly, we can conclude that Tdis ~ (0.05 0.1) J, but 
cannot provide a more accurate estimate. 

Note that the formation of a domain- wall network lead- 
ing to the intermixing of all six striped vortex patterns 
is impossible without the proliferation of free fractional 
vortices. This means that there is no possibility for the 
phase transitions related with the removal of the triple 
degeneracy of the uniaxial network states and suppres- 
sion of the algebraic phase behavior of Cg (r) to take place 
at different temperatures. On the other hand, the appear- 




FIG. 11: A boundary between two uniaxial networks with 
different preferable orientations of domain walls. 



ance of free fractional vortices leads to the screening of 
the logarithmic interaction of the conventional vortices, 
so that they also unbind. Accordingly, at T > Tdis, the 
system is in the completely disordered phase in which all 
correlations decay exponentially. 



VIII. CONCLUSION 
A. Phase diagram 

The main results of this work are summarized in Table 
I reviewing the properties of various phases represented 
in the phase diagram of the fully-frustrated XY model 
on a honeycomb lattice. The order in which these phases 
are listed corresponds to the increase in temperature. By 
the existence of the long-range order (LRO) in terms of 
vortex pattern we mean the selection of one of the six 
equivalent striped vortex patterns, see Fig. [T](a). On the 
other hand, the existence of the long-range order in terms 
of the domain- walls orientation implies the choice of one 
of the three preferable orientations of domain walls. The 
last three columns describe the behavior of the gauge- 
invariant phase correlation functions defined by Eq. ([7]) . 

The first line in Table I refers to the behavior of the 
correlation functions when the averages are calculated 
over the set of the all ground states of the model. A typ- 
ical ground state incorporates an irregular sequence of 
parallel straight domain walls, which leads to the expo- 
nential decay of the chirality correlation function, Eq. 
(|3]), as well as of phase correlation functions Cp(r) with 
p < 8 for all directions that are not exactly perpendicular 
to one of the directions of lattice bonds. Nonetheless, the 
system has the long-range order in terms of domain- walls 
orientation 

At T > domain walls acquire a finite free energy per 
unit length given by Eq. ([8|) and originating from the dif- 
ference in the free energy of spin waves.— At the lowest 
temperatures this leads to the suppression of infinite do- 
main walls crossing the whole system and the appearance 
of the long-range order in term of vortex pattern (phase 
No. 1 in Table I). The system spontaneously chooses one 
of the six equivalent striped vortex patterns. On the other 
hand, all phase correlation functions decay algebraically 
due to the presence of spin waves. 

With the increase in temperature, the first phase tran- 
sition takes place at T = Tdddw ~ 0.8 x 10~^ J and is 
related to the appearance of a sequence of fluctuating 
double-spaced double domain walls. On the average, the 
direction of these fluctuating walls is perpendicular to 
the direction of vortex stripes. This phase transition is 
continuous with the Ising-type behavior in a very nar- 
row critical region and the Pokrovsky-Talapov behavior 
in a wider region around the critical one. Although the 
transition taking place at T = Tdddw is related to the 
appearance of domain walls, above it (phase No. 2 in Ta- 
ble I) the discrete degeneracy (related to the long-range 
order in terms of vortex pattern) remains the same (six- 
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fold), whereas the decay of the phase correlation function 
Ci(r) changes from algebraic to exponential. This hap- 
pens because the presence of a doubly-spaced double do- 
main wall between two points shifts the phase difference 
by TT. 

The next phase which appears with the increase in 
temperature (phase No. 3 in Table I) has the same sym- 
metry properties as a typical ground state of the model. 
In particular, in this phase there is no long-range order in 
terms of vortex pattern. The main difference with a ran- 
dom sequence of parallel straight domain walls (charac- 
teristic for a typical ground state) consists in the presence 
of local defects which shift the positions of these walls 
but do not change their preferable orientation, see Fig. 
[Tjb). The phase transition to this uniaxial domain- wall 
network has to be of the first order and can be expected 
to take place when T is around 2rDDDW- The important 
feature of the uniaxial domain-wall network is that the 
nodes of this network, which in terms of phase variables 
correspond to logarithmically interacting fractional vor- 
tices, are bound in neutral pairs. 

The last phase transition leading to the complete sup- 
pression of the long-range order and to exponential de- 
cay of all correlation functions can be associated with the 
appearance of an isotropic domain- wall network in which 
some fractional vortices (the nodes of the network) arc 
not bound in pairs but free. By analyzing how the free 
energy of such a network depends on typical distance 
between free fractional vortices it turns out possible to 
establish that this transition must be of the first order 
and has to occur at Tdis ~ (0.05 ~ 0.1) J. 

All phase transitions mentioned above take place at 
temperatures which are much lower than the tempera- 
ture Tbkt ^ J at which the unbinding of pairs of the 
conventional vortices (with integer topological charges) 
would take place if any other fluctuations were absent. 
The existence at such temperatures of an exponentially 
small concentration of small bound pairs of integer vor- 
tices does not lead to any noticeable consequences. On 



the other hand, at T > Tdis the logarithmic interaction 
of integer vortices is screened due to the presence of a 
finite concentration of free fractional vortices and there- 
fore there are no grounds to expect the existence of an 
additional phase transition related to the unbinding of 
integer vortices. 



B. Discussion 

Thus we have established that the equilibrium ther- 
modynamics of the fully frustrated XY model on a hon- 
eycomb lattice supports the existence at T > of four 
different phases three of which are characterized by the 
presence of a long-range order. However the observation 
of some of these phases (and phase transition between 
them) in a real or numerical experiment may turn out to 
be a big problem. 

In particular, at < T < TdddWj the complete sup- 
pression of domain walls crossing the whole system (in- 
duced by the spin- wave contribution to their free energy, 
see Sec, imp is achieved only in the thermodynamic limit, 
^max — >■ oo. In a finite system with linear size imaxi in 
order to have the average number of such walls smaller 
than one, the size of the system should be larger than 
Lc(T), the temperature-dependent solution of the equa- 
tion 

imaxexp = 1 (18) 

At the point of the phase transition related to the ap- 
pearance of infinite double-spaced double domain walls, 
Lc ~ 3 X 10^. This means that numerical simulations of 
systems with L < 10^ does not allow one to observe 
the phase with suppressed domain walls independently 
of whether it is possible or not to reach thermal equilib- 
rium at T ^ Tdddw- 
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cxp 
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uniaxial domain-wall network 




+ 


exp 


exp 


a 


4 


disordered phase 






exp 


exp 


exp 



TABLE I: The properties of different phases of the fully frustrated XY model on a honeycomb lattice. In the last three columns 
"exp" denotes an exponential decay of the corresponding correlation function (for a generic direction) and "a" an algebraic 
decay. 



16 



More or less the same criterion applies to the observa- 
tion of the phase with a sequence of fluctuating double 
domain walls. Accordingly, for L < 10^ the increase of 
temperature will lead just to a smooth crossover from a 
random sequence of parallel straight domain walls to the 
uniaxial domain-wall network with the same preferable 
orientations of domain walls. However, a discontinuous 
transition from the uniaxial domain-wall network to an 
isotropic one (the disordered phase) can be expected to 
be observable even if the size of the system is not so 
large. On the other hand, it follows from the numerical 
simulations of Ref. [13 that the task of reaching thermal 
equilibrium at T ^ 0.1 J may require special efforts for 
its solution. 

In the case of the Berezinskii- Villain interaction the 
free-energy of spin waves is exactly the same for all 
domain-wall configurations and, therefore, does not lead 
to the removal of the accidental degeneracy of the ground 
states. In such a situation, the scenario of a smooth 
crossover from a random sequence of parallel straight do- 
main walls to the uniaxial domain-wall network works 
not only at imax ^ Lc bit also in the thermodynamic 
limit. That is, instead of three phase transitions an infi- 
nite system with the Berezinskii- Villain interaction must 
experience only one, related to the loss of the long-range 
order in the orientation of domain walls and of the alge- 
braic decay of C%{t). The conclusion that this transition 
has to be of the first order is not sensitive to a particular 
form of the interaction in the Hamiltonian of the model. 

Since the fully-frustrated XY model with the 
Berezinskii- Villain interaction on a honeycomb lattice is 
exactly equivalent to the half-integer Coulomb gas on a 
triangular lattice, the scenario described in the previous 
paragraph has to be realized also in such a Coulomb gas. 
Numerical simulations of the half-integer Coulomb gas on 
a triangular lattice were undertaken by Lee and Teitcl,^^ 
who have discovered a relatively sharp jump of the di- 
electric function at TjJ^ 0.04 but assumed it to be a 
finite-size induced artefact. This conclusion was based on 
the observation that the formation and motion of finite- 
energy local defects (analogous to those shown in our 
Fig. [5]) destroys the long-range order in vortex pattern at 
an arbitrary low temperature. However, this argument 
missed a possibility of having the long-range order in 
domain-walls orientation, which according to our anal- 
ysis should be present in the half-integer Coulomb gas 
on a triangular lattice at low enough temperatures. 

It should be mentioned that a uniformly frustrated XY 
model with this or that interaction provides just an ideal- 
ized description of a Josephson-j unction array or of a su- 
perconducting wire network. In physical situations there 
can exist other mechanisms for the removal of the ac- 
cidental degeneracy of the ground states related to the 
interactions not taken into account in the framework of 
an XY model. One of them is the magnetic interaction of 
currents in the junctions or in the wires. However, it is 
likely that the sign of the domain- wall energy induced by 
this interaction may depend on the particular geometry 



of the system. 

If -EdWi the domain wall energy per unit length, is 
positive, this will substantially improve the possibilities 
for the observation of the phase with the striped vortex 
pattern. Even if -Edw is few orders of magnitude smaller 
than J, it will be many orders of magnitude larger than 
the value of /dw given by Eq. ([5]) at the corresponding 
temperatures. The opposite sign of -Edw will lead to the 
stabilization of another periodic vortex pattern shown in 
Fig. 1(c) of Ref. [l3. In any case, a nonzero value of E'dw 
will allow the observation of a phase with a long-range 
order in terms of vortex pattern in the systems of less 
than macroscopic sizes. 

On the other hand, the magnetic interaction of cur- 
rents leads to the screening of the logarithmic interaction 
of vortices (both conventional and fractional) at large dis- 
tances. This will transform the phase transition related 
to the appearance of a sequence of fiuctuating double- 
spaced double domain walls into a crossover, because 
when the logarithmic interaction of fractional vortices is 
screened, such double walls no longer are topologically 
stable defects and can have free end-points. A substan- 
tial increase of i?DW may lead also to the disappearance 
of the region of stability of the uniaxial network state, 
which implies a direct transition from the phase with the 
long-range order in terms of vortex pattern into the dis- 
ordered phase. 

Another mechanism for the stabilization of striped vor- 
tex patterns in magnetically frustrated superconducting 
wire networks may be related with the non-uniformity of 
the order parameter amplitude^^ (see also Ref. jsj) . How- 
ever the conclusions of Ref. IBS' are in some contradiction 
with the results of Ref. [3^ devoted to including into anal- 
ysis the higher-order terms of the Ginzburg-Landau ex- 
pansion, which suggests that the problem requires further 
investigation. 

The uniformly frustrated XY model with triangular 
lattice and f = \ also allows for the formation of zero- 
energy domain walls parallel to each other,— so its phase 
diagram may have some common features with the one 
constructed in this work. In the uniformly frustrated XY 
model with / = ^ on a triangular lattice^^, as well as 
in the fully frustrated XY model on a dice lattice^ii^S 
the situation is more complex because these models al- 
low for the formation of two sets of parallel zero-energy 
domain walls which can cross each other without paying 
any energy for this. However, in such a situation it is also 
possible to expect the existence of a phase in which vor- 
tex pattern is disordered whereas the long-range order is 
related to the orientation of domain walls. 

On the experimental side, the conclusions of this work 
may be applicable not only to Josephson-junction arrays 
and superconducting wire networks with a half-integer 
number of flux quanta per plaquette, but also to mag- 
netically frustrated triangular arrays of microholeg'^'' or 
nanoholcs'^^ in thin superconducting films. For a half- 
integer number of fiux quanta per hole, the energy of 
different vortex configurations in such objects has to be 
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described by the Hamiltonian of the half-integer Coulomb 
gas with screened logarithmic interaction. 
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Appendix A: Berezinskii- Villain interaction and the 
Coulomb gas representation 

A magnetically frustrated network formed by identical 
superconducting wires can be described in London regime 
(when the amplitude of the superconducting order pa- 
rameter is uniform along the wires) by the Hamiltonian^ 



H 



J 



E 



|2 

Uj') 



(Ala) 



where djji is the integral of the gauge- invariant phase gra- 
dient along the wire connecting nodes j and j' . Variables 
Ojji can acquire arbitrary values, — oo < Ojjr < +00, but 
on any plaquette of the lattice have to satisfy the con- 
straint 



□ a 



03 



-27r/ (mod 2tt) , 



(Alb) 



where parameter / depends on the applied magnetic 
field and is equal to the number of superconducting fiux 
quanta per plaquette. 

In terms of phase variables ipj defined on the nodes 
of the network, the partition function of the model 
(IA1[) can be rewritten^ as the partition function of 
the uniformly frustrated XY model with the so-called 
Berezinskii- Villain interaction Vbv (^) defined by the re- 
lation 



exp 



Fbv(^) 



T 



E 



h—~oo 



(A2) 



The function Vbv(^) has the same symmetry and pe- 
riodicity as V{9) = —Jcos9. For J T it is every- 
where except the close vicinity of the point 9 = tt close 
to parabola, Vbv(^) ~ (J/2)0^ -I- const. In the opposite 
limit, J <^ T, the function Vbv(^) with an exponential 
accuracy is reduced to — Jeff cos 9 + const, where however 
coupling constant Jeff = 2Texp(— r/2J) is much smaller 
than J. 

The interaction function defined by Eq. (IA2[) was in- 
troduced by Berezinskii^ and Villain^ because it allows 
to simplify the analytical analysis of XY models. Namely, 
the integration over variables ipj in the partition function 
of an XY model with such an interaction is Gaussian and 
therefore can be performed exactly. In particular, in the 
case of the periodic boundary conditions the application 
of this procedure transforms the partition function of 



a uniformly frustrated XY model with the Berezinskii- 
Villain interaction into that of a "fractional" Coulomb 
gas^ described by the Hamiltonian 



HCG = ^'^ITlaGa 



(A3) 



where variables rua (the charges of the Coulomb gas) 
are defined on the sites a of the dual lattice and acquire 
values shifted with respect to integers by — /. Each of 
these variables is proportional to the sum of variables 
9jji --hjji (defined on lattice bonds) over the perimeter of 
plaquette a and, in terms of the XY representation, can 
be identified with the vorticity of this plaquette divided 
by 2n. 

The long-range interaction Gq^ entering Eq. (IA3p has 
a form 



(A4) 



where A is the discrete Laplace operator on the dual 
lattice. This interaction logarithmically depends on 
= Itq — r^l, the distance between a and /?, 



Go 



Gap = 27rro (Inr^^ -f k) 



(A5) 



where Fq is proportional to J but also depends on the 
structure of the lattice. In terms of the XY representation 
Fq is the helicity modulus describing the rigidity of the 
system with respect to phase twist. In the considered case 
To = J/VS, whereas k = ir/^/S for = 1 (the nearest 
neighbors on the triangular lattice) and monotonically 
increases by less than 0.5% with the increase of r^^ to 
infinity. 

It can be expected that the number of qualitative fea- 
tures of frustrated XY models with the conventional and 
with the Berezinskii- Villain interactions are the same. 
However, in situations when the ground states of the 
model possess an accidental degeneracy, one should be 
cautious because in a model with the Berezinskii- Villain 
interaction the free energy of the spin waves is exactly 
the same for all vortex configurations and therefore does 
not lead to the removal of such a degeneracy (in contrast 
to the models with the conventional form of the interac- 
tion) . 

In the case of a fully- frustrated model (with / = ^), 
the charges of the Coulomb gas are half-integer. Since 
the charges of the opposite signs attract each other it is 
rather clear that when the dual lattice is bipartite (square 
or honeycomb) the minimum of energy is achieved in 
configurations with a regular (checkerboard-like) alterna- 
tion of positive and negative charges. This is in a perfect 
agreement with having a two-fold discrete degeneracy of 
the ground states in terms of the XY representation. In 
the case of the half-integer Coulomb gas on a trangu- 
lar lattice it is impossible to construct a configuration in 
which each charge has neighbors only of the opposite sign. 
Accordingly, it is much less evident what are the states 
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minimizing the energy and how high is their degeneracy. In Appendix |B] we propose a transformation which sub- 
stantiahy simphfies this task and also provides a transparent expression for the energies of the excited states of the 
half-integer Coulomb gas. 



Appendix B: Effective charge representation for the energy of the half-integer Coulomb gas 

When considering the half-integer Coulomb gas described by Hamiltonian (jA3[) it is convenient to introduce auxiliary 
variables qa linearly related to ttIq, 

= + XAa/jinp , (Bl) 

where is the discrete Laplacian on the lattice on which variables ttIq are defined and the choice of A depends on 
the structure of this lattice. The replacement of rua by qa — X^apmp allows to rewrite Eq. (jA3l) as 

Hcg = \Y. lo'Gcfiqp + 47r2 J (A - A^z) ^ + 2^Va2 ^ (m„ + m^-f , (B2) 

Q,/3 a [aa') 



where z is the coordination number of the lattice, whereas 
the summation in the last term is performed over all pairs 
of nearest neighbors (aa'). The form of Eq. (IB2p suggests 
that the long-range interaction of the charges of the half- 
integer Coulomb gas can be replaced by their repulsion 
on neighboring sites plus the long-range interaction of the 
"effective charges" qa defined by Eq. (jBl[) . Naturally, it is 
convenient to choose the value of A in such a way that in 
the ground states all effective charges are equal to zero, 
which substantially simplifies the calculation of energy 
and also minimizes the first term in Eq. (jB2p . 

For A < 1/z the second term in Eq. (|B2p is minimized 
when nia = ±i, whereas the third term requires to max- 
imize the number of bonds on which m„ and vria' have 
opposite signs. On bipartite lattices it is always possible 
to do this for any pair of neighboring sites. In such a 
case, the fulfillment of condition = is achieved when 
one takes A = l/(2z). For square lattice this gives A = | 
and for honeycomb lattice A = i. For these values of A 
the states with the regular checkerboard alternation of 
positive and negative charges nia = ±1 ensure the si- 
multaneous minimization of all three terms in Eq. (jB2p 
which rigorously proves that they are the ground states. 

In the case of a triangular lattice, each plaquette has to 
have either one or three bonds with the same sign of nia 
and nia' ■ As a consequence of this, on the average, each 
charge has to have at least two neighbors of the same 
sign. In such a situation it is convenient to take A = ^, 
which for all configurations with = ±^ reduces Eq. 
(|B2|l to 

HcG ^\Y.^qo.G amp + [n + g , (B3) 

where 

(B4) 



N is the total number of sites and is the number of 
triangular plaquettes on which all three charges have the 
same signs. It is evident that the first term in Eq. (|B3[) is 
minimized when all variables qa are equal to zero and the 
second one when = 0. Since both these conditions can 
be fulfilled simultaneously, constant g is nothing else but 
the ground-state energy per site. Note that the value of 
g given by Eq. (jB4| is in perfect agreement with having 
in the fully frustrated XY model Ojj' = on one third of 
the bonds of the honeycomb lattice and Ojji = ±7r/4 on 
all other bonds. 

Thus, the set of the ground states of the half- integer 
Coulomb gas includes all states in which all charges are 
equal to ±1/2, each charge has exactly two neighbors 
of the same sign and on each triangular plaquette there 
are always two charges of one sign and one charge of the 
opposite sign. These states can be put into the correspon- 
dence with the ground states of the fully frustrated XY 
model on a honeycomb lattice by interpreting the charges 
nia as vorticities (divided by 27r) of the corresponding 
plaquettes, as it could be expected from the equivalence 
between the half-integer Coulomb gas and the fully frus- 
trated XY model with the Berezinskii- Villain interaction 
on the dual lattice (see Appendix However, the reduc- 
tion of Hamiltonian (jA3|) to the form (|B3P has allowed 
us to find the structure of the ground states of the half- 
integer Coulomb gas on a triangular lattice directly in 
terms of the Coulomb gas representation. 

Note that the form of Eqs. ([62]) and (|B3|l relies on 
the very special form of the interaction of charges in the 
Coulomb gas [given by Eq. (|A4I) ]. As soon as the form of 
the interaction is modified (for example, to take into ac- 
count the screening effects in a superconducting array or 
wire network), the possibility to transform the Hamilto- 
nian to the form (|B2[) with only local interaction of vari- 
ables ma disappears, which leads to a partial removal 
of the degeneracy between the states with qa ^ and 
A^3 =0. 
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In the configurations with nia = ±i on a triangular 
lattice, the values of effective charges Qa are given by 

^ ^.>, , (B5) 

where is the number of the nearest neighbors of site 
a with the same sign of the charge m. It follows from 
Eq. (jBSp that the effective charges of various sites are 
either zero or equal to ±i or multiples of these values. In 
Sec. IIVBI the appearance of logarithmically interacting 
vortices with fractional topological charges on those pla- 
quettes of the honeycomb lattice which have the number 
of neighbors with the same sign of vorticity not equal 
to two has been discussed directly in terms of the XY 
model. 

Eq. (jB3|) provides a convenient way for calculating the 
energies of local defects in half-integer Coulomb gas on a 
triangular lattice (or in the equivalent fully frustrated XY 
model with the Berezinskii- Villain interaction on a hon- 
eycomb lattice) by summing just few terms. The values of 



these energies for various finite-energy defects discussed 
above and illustrated in Figs. are listed in Table II. 
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TABLE II: The energies (in units of J) of the finite-energy 
defects presented in Figs. [3][9] for the case of the Berezinskii- 
Villain interaction. 



S. Teitel and C. Jayaprakash, Phys. Rev. Lett. 51, 1999 
(1983). 

M. Hasenbusch, A. Pelissetto and E. Vicari, J. Stat. Mech.: 



Th. Exp. P12002 (2005) , [ajOCivxond-ma t /0509682 (2005). 
S.E. Korshunov, Usp. Fiz. Nauk 176, 233 (2006) [Phys. 
Usp., 49, 225 (2006)]. 

For reviews, see R. S. Newrock, C.J. Lobb, U. 
Geigenmiiller and M. Octavio in Solid State Physics, 
Vol. 54, edited by H. Ehrenreich and F. Spaepen (Aca- 
demic Press, San Diego, 2000), p. 263 and P. Martinoh 
and Ch. Leemann, J. Low Temp. Phys. 118, 699 (2000). 
J. Villain, J. Phys. C 10, 1717 (1977). 
L. B. loflfe and M. V. Feigel'man, Phys. Rev. B 66, 224503 
(2002); B. DouQot, M.V. Feigel'man and L.B. loflfe, Phys. 
Rev. Lett. 90, 107003 (2003). 

V. L. Berezinskii, Low-temperature properties of the two- 
dimensional systems, thesis (L. D. Landau Institute, 1971), 
pubhshed by Fizmatlit, Moscow, 2007 (in Russian). 
J. ViUain, J. Phys. (Paris) 36, 581 (1975). 
E. Fradkin, B. Huberman and S. H. Shenker, Phys. Rev. 
B 18, 4789 (1978). 

S.E. Korshunov, Phys. Rev. Lett. 88, 167007 (2002). 

P. Olsson, Phys. Rev. B 55, 3585 (1997). 

W.Y. Shih and D. Stroud, Phys. Rev. B 32, 158 (1985). 

S.E. Korshunov, J. Stat. Phys. 43, 17 (1986). 

S.E. Korshunov and B. Dougot, Phys. Rev. Lett. 93, 

097003 (2004). 

J.-R. Lee and S. Teitel, Phys. Rev. B 46, 3247 (1992). 
W.Y. Shih and D. Stroud, Phys. Rev. B 30, 6774 (1984). 
R.W. Reid, S.K. Bose and B. Mitrovic, J. Phys.: Cond. 
Matt. 9, 7141 (1997). 

R. N. Bhatt and A. P. Young, Phys. Rev. B 37, 5606 (1988). 
J. Villain, R. Bidaux, J. P. Carton and R. Conte, J. Phys. 
(France) 41, 1263 (1980). 

E.F. Shender, Zh. Eksp. Teor. Fiz. 83, 326 (1982) [Sov. 
Phys. - JETP 56, 178 (1982)]. 

H. Kawamura, J. Phys. Soc. Jpn. 53, 2452 (1984); S.E. 



Korshunov, Pis'ma Zh. Eksp. Teor. Fiz. 41, 525 (1985) 
[JETP Lett. 41, 641 (1985)]; J. Phys. C 19, 5927 (1986); 
C.L. Henley, J. Appl. Phys. 61, 3962 (1987); Phys. Rev. 
Lett. 62, 2056 (1989). 

S.E. Korshunov, A. Vallat and H. Beck, Phys. Rev. B 51, 
3071 (1995). 

Analogous hidden symmetry exists also in the fully frus- 
trated XY model on a dice lattice, see Ref. 24. 
S.E. Korshunov, Phys. Rev. B 71, 174501 (2005). 
B. D. Metcalf, Phys. Lett. 46A, 325 (1974); M. Kaburagi 
and J. Kanamori, Jpn. J. Appl. Phys., Suppl.: pt. 2, 145 
(1974). 

S.E. Korshunov, Phys. Rev. B 72, 144417 (2005). 
^'^ S. E. Korshunov, Pis'ma Zh. Eksp. Teor. Fiz. 41 216 (1985) 
[Sov. Phys. - JETP Lett., 41, 264 (1985)]; J. Phys. C 19, 
4427 (1986); D. H. Lee and G. Grinstein, Phys. Rev. Lett. 
55, 541 (1985). 

V. L. Pokrovsky and A. L. Talapov, Zh. Eksp. Teor. Fiz. 
78, 269 (1980) [Sov. Phys. - JETP, 51, 134. (1980)]; V. L. 
Pokrovsky and A. L. Talapov, Theory of Incommensurate 
Crystals (Harwood Academic, Chur, 1984). 
T. Bohr, V. L. Pokrovsky and A. L. Talapov, Pis'ma Zh. 
Eksp. Teor. Fiz. 35, 165 (1982) [Sov. Phys. - JETP Lett., 
35, 203 (1982)]; T. Bohr, Phys. Rev. B 25, 6981 (1982); 
H.J. Schulz, Phys. Rev. B 28, 2746 (1983). 
When the helicity modulus is an anisotropic tensor, F in 
Eq. p6p stands for the geometric mean of its eigenvalues. 
This estimate does not take into account the renormal- 
ization effects (related to the presence of bound pairs of 
fractional vortices), which would shift Tfv further down. 
J. Villain, in Order in Strongly Fluctuating Condensed 
Matter Systems, edited by T. Riste (Plenum, New York - 
London, 1980), p. 221; J. Villain, Surf. Sci. 97, 219 (1980); 
J. Villain and M.B. Gordon, Surf. Sci. 125, 1 (1983). 

^3 Y. Xiao, D. A. Huse, P. M. Chaikin, M. J. Higgins, S. Bhat- 
tacharya and D. Spencer, Phys. Rev. B 65, 214503 (2002). 

^'^ O. Sato and M. Kato, Physica B 329-333 1405 (2003). 



20 



P. Erdos and W. Zheng, Phys. Rev. B 82, 134532 (2010). 

S.E. Korshunov, Phys. Rev. B 63, 134503 (2001). 

T. Ishida, M. Yoshida, K. Okuda, S. Okayasu, M. Sasasc, 

K. Hojou, A. Odawara, A. Nagata, T. Morooka, S. 

Nakayama and K. Chinone, Physica C 357-360, 604 

(2001). 

M. D. Stewart, Jr., Z. Long, J.M. Valles, Jr., A. Yin and 
J.M. Xu, Phys. Rev. B 73, 092509 (2006); H. Q. Nguyen, 



S. M. HoUcn, M. D. Stewart, Jr., J. Shainlinc, A. Yin, J. M. 
Xu, and J.M. Valles, Jr., Phys. Rev. Lett. 103, 157001 
(2009). 

S. Alexander, J. Physique 44, 805 (1983). 

A. Vallat, S. E. Korshunov and H. Beck, Phys. Rev. B 43, 

8482 (1991); A. Vallat and H. Beck, Phys. Rev. B 50, 4015 

(1994). 



